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ABSTRACT 


The theory of the solar gravitational figure is derived including the 

effects of differential rotation. It is shown that is smaller than J2 by a 

2 

factor of about 10 rather than being of order (J2) as would be expected for 
rigid rotation. The dependence of both J2 and on envelope mass is given. 

High order p-mode oscillation frequencies provide a constraint on solar structure 
which limits the range in envelope mass to the range 0.01 < M_/M < 0.04. For 

C V 

an assumed rotation law in which the surface pattern of differential rotation 
extends uniformly throughout the conyCctive envelope, this structural constraint 
limits the ranges of J2 and in units of 10 

< 1.5. Deviations from these ranges would imply that the rotation law is not 
constant with depth and would provide a measure of this rotation law. 
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"This work was performed for the Jet Propulsion Laboratory, California 
Institute of Technology sponsored by the National Aeronautics and 
Space Administration under Contract NAS7-100." 
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to 10 < J2 < 15 and 0.6 < - 
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TIIK SOLAR GRAVITATIONAL FIGURE — AND 
R.K. Ulrich and G.W. Hawkins 


I. INTRODUCTION 


Tlie gravitational Figure oF the sun can be measured directly by sending 
a probe carrying a drag-Free guidance system to within a Few solar radii of the 
solar center. The interpretation of such a measurement will require knowledge 
of how the solar structure and rotational law influence the multi-pole expan- 
sion of the external gravitational field. The study of solar oscillations limits 
the variations in due to structural uncertainties to less than a factor of 2. 
Refinements in the observed frequencies of the non-radial modes of high order can 

further reduce this uncertainty. In addition, v;e show how the existence of dif- 

2 

ferential rotation throughout the convective zone will produce >> An 

accurate measurement of J2 and detection of can substantially increase our 
knowledge of the solar rotation law. 

The papers dealing with the theory of the gravitational figure for rotating 
stars dlyide into two types. ‘ The theory of rapidly rotation configurations is 
generally limited to the study of polytropic models. The classical Clalrant- 
Legendre expansion resulting from this approach is well-describe by Tassoul 
(1978) . Applications based on this metliod are discussed by James (196A) , 

Ostriker and Mark (1968) and Hubbard, Slattery and DeVito (1975). While this 
method is very powerful in its ability to treat models of objecits in wlilch con- 

(I \ ' 

trlpetal acceleration and gravitational acceleration approach equality, it is 
not easily applied to non-poly.troplc modelii because of the fundamental way in 
which the poly tropic assumption is used. 
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I’he suGoiul type of approach Is that described by Goldrelch and Schubert 
(1972). In this approacli* the stellar structure equations are analyzed in a 
perturbation expansion. Consequently, this method fails for objects rotating 
fast enough tltat the centripetal and gravitational accelerations are comparable. 

It lias tlie advantage for application to tine solar case that it is not restricted 

\\ 

to poly tropic models. Since the sun is rotating slowly, we use the Goldreicii 
and Schubert (1972) approach. We present in the following section tlie theoreti- 
cal analysis of a dlffereuijl.aliy rptating model and derive the equation which 
doscrihes the octopnle moment of the sun, In i III we describe the nop- 

rotating zero order models and present the results for J 2 and in § XV. 

II. THEORY OK J, 

We adopt os a simple rotation law: 

n = (1) 

wliere 0^^ is the equatorial angular rotation rate and p is the cosine of tlie 

2 

colatitude. The term f2,> p represents the differential rotation. Fitting the 

rotation law of Howard, Boyden and Laboute (1980) to Equation (1), we find that 

~ “0*7. We assume that ^2 constant in tha/^onvective envelope 

and that = 0 below the convective envelope. Although tli^s latter assumption 

•may not be valid, 02 n'Ust go to zero at r = 0 because otherwise there would be 

2 

a divergence in tlie shear velocity field. The centripetal acceleration r 0 

2 2 2 

appears in tlie equations of momentum balance. We expand 0 as 0^ (1 + a p ) 

(( ^ ^ 

vvhe^re: * « 
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Following Goldreich and Schubert (1972), we write the momentum balance equa- 
tion in spherical coordinates (r, p) in terms of the gravitational potential <[',• 
the pressure? p and the density p. The r and p com|)onents of tills equation aro: 


0 IJ- " “ P Oq (1 + ap^) r (1 - p^) 

p ^ a 4. p jj2 ^,2 ^ g^2j ^ 
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Assuming that the acceleration of the rotation is small compared to the accelera- 
tion of gravity, we may expand all quantities into unperturbed quantities (^Q, p^ 
and p^ and tlieir perturbations <|i|> Pj^ end pj^. The unperturbed equations are stan- 
dard and will not be repeated here. The first-order perturbation eq^tions are: 
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We now expand p^^ and Pj^ in spherical harmonics, e.g.: 


+ p„ r^ (1 + ap^) p. 
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h ” A^O.2,4 


(5) 

( 6 ) 


(7) 


where P^(p) Js the Legendre polynomial of order £. Using standard recursion 
relations for the P , the momentum equations have S, = 2 and Z 1 = 4 components 

• ■ !i ■ 

which are: / 
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The above equations must be supplemented by tlie Poisson equations 
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Equations (8) and (9) and the Jl ® 2 component of Equation (12) combine to 
give the analogue of tlie standard. Goldreich and Schubert (1972) equation for 
This equation iss 
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This differs from the standard equation in having the extra term involving 4/21 


and in replacing Oq with Oq( 1 + 3a/7)^ The equation for the 2=4 component of 


Equation (12) is similarly derived and iss 
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iSqua cions (13) and (1^*) must be solved subject to the boundary conditions 
j) — Jl-1 

chat at r " 0 and “ r as r « Numerically, theue conditions 

j_ 

were satisfied by expanding as r |it the innermost mass shell. Equations 
(13) and (14) were then djctegrat from ^ihe center to the outermost mass point. 

The outer boundary condition was imposed as an equation of condition and was 

adjusted through a Newton-Raphson iteration procedure. The distribution of a(r) 

was chosen so that a(r) 0 for r < r and a(r) « - 1.4 for r > r where r is 

C C (3 

the radius at Che inner edge of the convection zone. In crossing this radius, 
the following jump condition is applied to the solution: 


!!ia I I ”0 "o " 

8r I r +c “ 9r • r -e ” 35 M„ 

e e r 


(15) 


where e is a small distance encompassing the gradient in a. The exterior gravi- 
tational field is described by the normalized values of 



i+1 


GM nir 


( 16 ) 


m. , THE SOLAR models 

- The static solar model has beik calculated using the input physics 
described in the paper by IJahcail etyal. (1980). The method differs front 
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pcdvioua solar model calculations that have been done usXnB the UCLA code In 
two respects! ' 

1) The atmosphere code discussed by Ulrich and Rltodes (1977), Ulrich, 

t 

Rhodes and Deubner (1979) and Lubow, Rhodes and Ulrich (19R0), has been 
extended to include all of the interior physics, dy adjusting the mix- 
ing length and hydrogen abundance, the atmosphere integration has been 
extended inward to a radius of 1,5 x 10^'^ cm, The final interior tone 

is treated analytically with a theory accurate to better than l/percent. 

/) ^ ^ [I ' 

■ Tlie interior distribution of hydrogen abundance has been taken Ifrom the 

current standard solar model. This method yields neutrino fluxes whicli 
exceed those in the standard model by 10 percent. 

2) Partial ionization of all heavy elements Is now included in an approxi- 
mate way following Bodenlielmer, Forbes, Gould and Henyey (1965). Fur- 
ther details of the calculational methods will be, given by Ulrlcli and 
Rliodes (1981). 

Two fully converged models were computed: Tlie standard model and a modal 

with the interior heavy element abundance reduced from Z <= 0.018 to Z » 0.005. 

The convective envelope abundances were normal. The mix of both cases was that 

given by Ross and Aller (1976). The second model is referred to as the low Z 

‘ I 

model. In addition to these models, a series of incomplete envelopes was ca^cu- 

- ' 

iatod with a specified mixing Icngtli parameter, ifnose envelopes were stoi/ped 
wlien either 99 percent of the mass or radius was used in tlie inward Integration 
leaving B large residual of the other variable. An approximate analytic continu- 
ation to the center was added by assuming that ■* r , (Pq - P) « r , (p^ - p) 


2 2 

« V and (Tq “ T) « r . Almost none of the perturbation in is produced In 
the poorly treated portion of the model. 

IV. RESULTS and DISCUSSION 

There are two factors which can Influence the values of 

distribution of mass through the sun and 2) the distribution of n (r, 0) through 

the solar interior. This latter function is observed at r ** R , but tliete is 

0 

as yet only a preliminary indication of its variation with r. A range in 
structure variation is possible in principle so that both factors must be con- 
sidered in interpreting any potential measurement of and J^. Fortunately, 
the structural variations are not constrained by the measurements of tlie eigen- 
frequencies of high order non-radlal oscillations. The recent observations by 
Rhodes, Harvey and Duvall (1981) limit the variations in the env^^lope mass to 
0.01 < M /M < 0.04. As a measure of the importance of structure in determining 

GO 

J 2 and we have assumed that D (r, 0) = D (R, 0) for r > r^ and that fj (r, 0) 
0 (R) for r < r^. We have used the rotation law given by Howard, Boydon and 
La Bonte (1980) in calculating J., and J^. Our results are given in; Figure 1. 

The above l^imits on translate into permitted ranges in J 2 and J^: 

1.0 X lO"^ < J, < 1.5 X 10 "^ 

.9 -8 

6 X 10 ^ < 1.5 X 10 

Heasured values of J 2 or outside of these ranges would indicate a deviation 
from our simple assumptions for 0 (r, 0). 

w 

The subsurface rotation rate can be measured, with the method described by 
Rhodes, Deubner and Ulrich (1979) and by Deubner, Rhodes and Ulrich (1979). 


This second paper gives some first results which show fi (r, 0) increasing as r 
decreases. If ((this result is borne out by subsequent observations! we would 
expect larger values for and The theory of the solar dynamo depends 

on such a variation of n with r (Stix 1977). Ultimately, It may be possible 
to probe the rotation law using the oscillations; however, there are two criti- 
cal uncertainties which must be resolved before we will know the effectlvness 
of the oscillation method*. 1) the oscillations must be adequately long lived 
and 2) it must be possible to identify correctly individual modes of oscill.n- 
tlon for all spatial scales. Also, the oscillations cannot provide a good 
measure of the rotation rate off the solar equator because foreshortening limits 
the viewing region. Thus, the measurement of both J 2 and would provide very 
important constraints on 0 (r, 0), which are different from the constraints 
imposed by the oscillations. 

We would like to thank Jim Underwood for Ills helpful conunents on Che manu- 
script. 
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FIGURIi CAPTION 


IMR ura 1 . The depenclence of J2 and on envelope mass 

Coprespondlng vftJuea of IfW ace indicated f»m the cop of the figure. 
Tliese values are dopondcnu on a numbor of coinputntlonnl decnlls and are 1 u«h 
significant than M„/hl . . The chnracCerlBticB of Che standard model and the low 
7, model are discussed In the text. 




